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HIGHLIGHTS 


•  A  formal  design  optimization  methodology  for  functionally  graded  electrodes  of  Li+  batteries  is  developed. 

•  A  multi-scale  finite  element  model  for  predicting  the  power-storage  characteristics  of  a  full  battery  cell  is  presented. 

•  A  multi-objective  formulation  is  introduced  to  optimize  the  useable  capacity  over  a  range  of  discharge  currents. 

•  The  influence  of  spatially  varying  porosity  and  sizes  of  active  particles  on  optimized  electrodes  is  analyzed. 

•  The  difference  between  considering  a  half-cell  versus  a  full  cell  in  the  optimization  process  is  shown. 
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Design  optimization  for  functionally  graded  battery  electrodes  is  shown  to  improve  the  usable  energy 
capacity  of  Li  batteries  predicted  by  computational  simulations  and  numerically  optimizing  the  electrode 
porosities  and  particle  radii.  A  multi-scale  battery  model  which  accounts  for  nonlinear  transient  trans¬ 
port  processes,  electrochemical  reactions,  and  mechanical  deformations  is  used  to  predict  the  usable 
energy  storage  capacity  of  the  battery  over  a  range  of  discharge  rates.  A  multi-objective  formulation  of 
the  design  problem  is  introduced  to  maximize  the  usable  capacity  over  a  range  of  discharge  rates  while 
limiting  the  mechanical  stresses.  The  optimization  problem  is  solved  via  a  gradient  based  optimization.  A 
LiMn204  cathode  is  simulated  with  a  PE0-LiCF3S03  electrolyte  and  both  a  Li  Foil  (half  cell)  and  LiC6 
anode.  Studies  were  performed  on  both  half  and  full  cell  configurations  resulting  in  distinctly  different 
optimal  electrode  designs.  The  numerical  results  show  that  the  highest  rate  discharge  drives  the  sim¬ 
ulations  and  the  optimal  designs  are  dominated  by  Li+  transport  rates.  The  results  also  suggest  that 
spatially  varying  electrode  porosities  and  active  particle  sizes  provides  an  efficient  approach  to  improve 
the  power-to-energy  density  of  Li+  batteries.  For  the  half  cell  configuration,  the  optimal  design  improves 
the  discharge  capacity  by  29%  while  for  the  full  cell  the  discharge  capacity  was  improved  61%  relative  to 
an  initial  design  with  a  uniform  electrode  structure.  Most  of  the  improvement  in  capacity  was  due  to  the 
spatially  varying  porosity,  with  up  to  5%  of  the  gains  attributed  to  the  particle  radii  design  variables. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Further  advances  in  the  energy  density,  rate  capability,  and 
cycling  efficiency  of  batteries  are  needed  in  order  to  meet  growing 
energy  storage  demands.  The  essential  characteristics  of  a  battery 
are  typically  represented  via  (1)  Ragone  plots  showing  the  tradeoff 
between  energy  storage  capacity  and  power  (or  rate)  capability  of 
the  battery  and  (2)  cycling  efficiency  values  which  measure  the  loss 
in  energy  storage  capacity  as  a  function  of  cycle  number.  In  this 
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paper,  we  focus  on  the  trade-off  between  power  and  energy  den¬ 
sities  with  consideration  of  mechanical  stress  levels  that  could  lead 
to  degradation  over  time.  We  do  not  explicitly  model  degradation 
phenomena  such  as  loss  of  lithium  due  to  side  reactions  and  SEI 
layer  formation  and  mechanical  fracture  [1];  instead  we  consider 
approximately  mechanical  degradation  by  introducing  a  stress 
limit  in  the  electrode  particles. 

In  general  the  energy  storage  capacity,  or  usable  capacity  de¬ 
creases,  with  increasing  rate  or  power  demand  as  seen  in  Fig.  1. 
Here,  we  define  the  usable  capacity  as  the  amount  of  energy  that 
can  be  extracted  from  the  battery  during  discharge.  In  this  study, 
we  will  show  that  the  decrease  in  usable  capacity  at  higher 
discharge  rates  can  be  mitigated  by  optimally  designing  the  inter¬ 
nal  electrode  structure  based  on  the  range  of  rates  required  from 


240 


S.  Golmon  et  al.  /  Journal  of  Power  Sources  253  (2014)  239-250 


Fig.  1.  Ragone  plot  comparing  battery  technologies  [34]. 


the  battery  for  a  given  application.  We  focus  the  study  on  lithium 
batteries  which  are  among  the  most  promising  battery  chemistries 
due  to  their  high  energy  densities  and  operating  voltages.  Design 
variables  for  the  battery  electrode  include  the  local  porosity  and 
particle  radii  distributions  in  the  electrodes.  Additional  design 
variables  such  as  the  electrode  thickness  and  material  properties 
can  be  incorporated  into  this  approach.  We  formulate  this  design 
problem  via  a  multi-objective  formulation,  maximizing  the  usable 
capacity  at  various  discharge  rates  while  constraining  the  stress 
levels  in  the  electrode  particle  to  be  below  a  specified  maximum 
value.  The  optimization  problem  is  solved  via  a  gradient-based  al¬ 
gorithm  which  allows  for  a  large  number  of  both  design  constraints 
and  optimization  variables.  The  battery  simulation  uses  a  transient, 
non-linear,  multi-scale  framework  to  simulate  the  discharge  of  the 
battery  at  prescribed  rates.  The  sensitivities  of  the  objective  and 
constraints  to  the  design  variables  are  calculated  via  a  multi-scale 
adjoint  approach  which  is  more  accurate  and  efficient  than 
computation  of  the  sensitivities  via  finite  differencing  [2].  Adjoint 
sensitivity  analysis  combined  with  a  gradient-based  optimization 
algorithm  drives  the  problem  to  the  optimal  electrode  structure. 

This  paper  is  aimed  at  advancing  the  state-of-the-art  in  model- 
based  optimization  for  Li  batteries.  Previous  optimization  studies 
for  Li  battery  performance  have  focused  on  pack-level  thermal  [3— 
5]  and  power  management  [6]  strategies.  Design  studies  focusing 
on  the  cell  level  have  been  severely  limited  in  number  of  design 
variables  due  to  the  complexity  and  computational  cost  of  most 
battery  models  [7-10].  This  brief  literature  review  does  not  include 
design  studies  based  on  experimental  investigations  and  parameter 
sweeps.  Here,  we  are  not  limited  in  number  of  design  variables  due 
to  our  use  of  multi-scale  adjoint  sensitivity  analysis  developed  by 
Golmon  et  al.  [2].  Our  approach  is  flexible,  allowing  for  different 
types  of  objectives,  constraints,  and  types  of  design  variables 
without  compromising  the  accuracy  of  the  overall  model.  Possible 
optimization  variables  include  but  are  not  limited  to:  the  electrode 
porosity  and  particle  size,  discharge  rate(s),  electrode  thickness, 
and  material  properties.  In  this  study,  we  focus  on  improving  the 
discharge  performance  of  the  battery  through  functionally  graded 
electrodes  by  optimizing  the  porosity  and  radii  distribution  in  both 
electrodes  over  a  range  of  discharge  rates.  We  consider  the  effects 
of  Li+  transport  in  the  electrolyte,  Li  transport  in  the  electrode 
particles,  local  electrode  surface  kinetics,  diffusion-induced 
stresses  in  the  electrode  particles  and  resulting  aggregate 
swelling  of  the  electrode.  Solid  electrolyte  interphase  (SEI)  layer 
growth,  chemical  side  reactions,  thermal  effects,  and  degradation 
effects  from  cycling  and  aging  are  not  included.  The  optimization 
approach  used  here  could  be  extended  to  more  sophisticated 
models  such  as  that  of  Kim  et  al.  [11]. 


2.  Multi-scale  battery  model 

A  typical  Li  battery  is  shown  in  Fig.  2  and  consists  of  current 
collectors,  porous  electrodes  and  a  separator  region.  The  porous 
electrode  includes  active  insertion  compounds,  conductive  addi¬ 
tives,  and  binders  with  electrolyte  filling  the  pores.  The  liquid 
volume  fraction  of  electrolyte  is  referred  to  as  the  porosity,  e,  and  in 
a  typical  battery  is  30%  [12].  When  the  battery  is  discharged,  Li  in 
the  anode  oxidizes.  Electrons  flow  through  the  external  circuit  to 
the  cathode  and  lithium  ions,  Li+  are  carried  through  the  electrolyte 
via  diffusion  and  migration  to  the  cathode,  where  they  are  reduced. 
In  modern  Li  batteries,  Li  inserts  into  a  host  electrode  material 
rather  than  plating  onto  metals,  and  these  compounds  shrink  and 
swell  as  a  result  of  the  Li  insertion  and  de-insertion.  The  resulting 
stresses  can  exceed  the  fracture  toughness  of  the  active  particle  and 
experimental  studies  have  linked  mechanical  effects  with  capacity 
fade  13]  and  localized  degradation  of  the  electrode  [14].  Stress- 
induced  surface  cracks  in  electrode  particles  have  been  shown  to 
nucleate  new  SEI  layer  growth  and  substantially  contribute  to  the 
loss  of  Li,  causing  capacity  fade  1  ].  These  studies  demonstrate  that 
it  is  necessary  to  consider  mechanical  effects  when  optimizing  the 
electrode  structure  in  order  to  limit  mechanical  degradation.  In  this 
study,  we  define  performance  in  terms  of  the  trade-off  between 
usable  energy  and  power  density  during  discharge  of  the  battery.  To 
minimize  mechanical  degradation  over  the  life  cycle  of  the  battery, 
we  limit  electrode  stress  levels  during  discharge. 

Here,  we  use  a  coupled  electrochemical-mechanical  multi¬ 
scale  model  to  simulate  the  discharge  of  Li  batteries  [15].  This 
model  is  an  extension  of  Doyle  and  Newman’s  electrochemical 
battery  model  [16-20]  which  describes  the  transport  processes 
within  the  electrodes  using  porous  electrode  theory  [21,20]  rather 
than  modeling  the  detailed  geometry  of  the  composite  electrode. 
We  have  extended  this  model  to  include  deformations  due  to 
external  mechanical  loads  and  swelling  of  electrode  particles  upon 
lithiation.  This  model  explicitly  considers  two  length  scales:  at  the 
macro-scale,  transport  processes  in  the  electrolyte,  electric  poten¬ 
tials  in  both  the  electrolyte  and  electrode,  and  mechanical  de¬ 
formations  across  the  battery  cell  are  described;  at  the  micro-scale, 
the  response  of  a  single  electrode  particle  is  modeled;  a  meso-scale 
model  is  used  to  relate  these  two  scales  using  volume  averaging 
homogenization  methods.  The  effective  macro-scale  properties  are 
computed  as  a  function  of  the  micro-scale  state  variables  at  each 
time-step  in  the  transient  simulation.  Qualitative  comparisons  of 
discharge  curves  generated  by  the  electrochemical  portion  of  this 
model  with  experimental  studies  have  shown  good  agreement 
[22]. 

In  the  following  sub-sections,  we  summarize  the  governing 
equations  used  to  describe  the  battery  behavior  at  the  macro-scale, 
micro-scale  and  the  homogenization  approaches  used  to  connect 
the  two  scales  at  the  meso-scale. 

2.1.  Governing  equations  at  macro-scale 

At  the  macro-scale,  the  Li+  concentration  in  the  electrolyte,  q, 
electric  potentials  in  the  solid,  and  liquid,  </>2*  phases,  and  dis¬ 
placements,  u,  are  modeled. 


Anode  Separator  Cathode 


Fig.  2.  Battery  cell. 


S.  Golmon  et  al.  /  Journal  of  Power  Sources  253  (2014)  239-250 


241 


Lithium  ions  are  transported  through  the  electrolyte  via  diffu¬ 
sion  and  migration: 


9C1  V, 

'8T+7 


N  +  F937i2-VC'-(1-t°)Jeff  =  °’ 


(1) 


where  N  is  the  diffusive  flux,  is  the  Li+  transference  number,  F  is 
the  Faraday  constant,  and  jeff  is  a  source  term  which  accounts  for 
the  volume  average  flux  of  Li+  from  the  electrolyte  into  the  elec¬ 
trode  particles. 

Within  the  solid  (electrode),  the  electric  potential  is  governed  by 
Ohm’s  law: 


Table  1 

Boundary  conditions  for  macro-scale  battery  model. 


Anode-current 
collector  boundary, 

r  ac 

Cathode-current 
collector  boundary, 
Fee 

Anode-separator 
boundary,  TA s 

Li+  flux 

V  Q  =  0 

V  Ci  =  0 

n  =  L1/1 

Current  in  solid 

it-/ 

ii  =/ 

ii  =  0 

particles 

Current  in 

i2  =  0 

i2  =  0 

i  2  =  1 

electrolyte 

Displacements 

u  =  0 

u  =  0 

u  =  0 

2.2.  Governing  equations  at  micro-scale 


V-ii  +Fjeff  =  0.  (2) 

Similarly,  the  electric  potential  in  the  electrolyte  is  governed  by 
a  modified  Ohm’s  law  that  accounts  for  the  dependence  of  the 
current  in  the  electrolyte,  i2,  on  the  concentration  of  Li+  in  the 
electrolyte,  c\. 

V-i2-Fjeff  =  0.  (3) 

The  macro-scale  mechanical  response  is  assumed  to  be  linear- 
elastic  with  a  volume-averaged  eigenstrain  model  to  account  for 
swelling  of  the  micro-scale  electrode  particles: 

V<r  +  b  =  0,  (4) 


At  the  micro-scale,  the  response  of  a  single  electrode  particle 
embedded  in  the  host  electrolyte  is  modeled.  We  assume  that  the 
macro-scale  Li+  concentration,  electric  potentials,  and  displacement 
fields  are  spatially  constant  in  the  electrolyte  immediately  sur¬ 
rounding  the  particle.  Assuming  the  micro-scale  particles  are 
spherical  results  in  uniform  boundary  conditions  and  allows  for 
simplification  of  the  micro-scale  model  to  one  dimension  using 
spherical  coordinates.  The  micro-scale  model  is  based  on  the  coupled 
diffusion-stress  model  developed  by  Zhang  et  al.  [23],  and  extended 
to  include  Butler— Volmer  surface  kinetics  and  surface  pressures. 

The  concentration,  cs,  and  deformation,  ur,  fields  within  the  par¬ 
ticle  are  described  by  diffusion  and  linear  static  mechanical  models: 


where  a  and  b  are  the  macro-scale  stress  tensor  and  body  force. 
Here,  the  body  force  is  assumed  to  be  zero. 

The  diffusive  flux,  electric  currents,  stress  and  strains  are 
described  by  the  following  constitutive  equations: 


N  =  — DeffVC|, 


(5) 


9cs 
"a t 


a 

ar 


J  =  0, 


-jr  +  j  (ffR  -  °t)  =  0, 


with  the  following  constitutive  equations: 


(9) 

(10) 


it  =  -AV0!, 


(6)  „  /  Qcs 

J  —  —Ds  (  VrCs  —  Vr(7h 


(ID 


h  = 


~Keff 


^(l-t°)vln(q) 


(7) 


<T  =  C  :  (e-ech),e  =  i(vu  +  Vur),  (8) 

where  Deff  is  the  effective  diffusion  coefficient  of  Li+  in  the  elec¬ 
trolyte,  A  is  the  conductivity  of  the  electrode  particles,  /ceff  is  the 
effective  ionic  conductivity  of  the  electrolyte,  Ceff  is  the  homoge¬ 
nized  elasticity  tensor,  and  eCh  is  the  volume-averaged  chemically 
induced  eigenstrain  due  to  the  aggregate  swelling  of  the  electrode 
particles.  More  complex  mechanical  models  to  capture  nonlinear 
stress— strain  relationships,  finite  strains,  and  irreversible  de¬ 
formations  could  be  included  in  this  framework. 

For  a  galvanostatic  discharge  of  the  battery,  boundary  conditions 
are  applied  as  follows:  at  the  anode-current  collector  boundary, 
Tac,  and  cathode-current  collector  boundary,  Tec,  all  the  current  is 
carried  by  the  electrode  particles.  The  boundary  conditions  are 
given  in  columns  two  and  three  of  Table  1. 

To  simplify  the  computational  burden  of  these  simulations,  we 
also  model  a  Li  half  cell,  in  which  a  porous  cathode  is  discharged 
against  a  Li  foil  anode;  the  anode  is  modeled  through  boundary 
conditions  at  the  anode-separator  interface.  At  the  anode-separator 
boundary,  Tas.  Li+  enters  the  electrolyte  based  on  the  rate  of 
discharge,  J,  and  all  current  is  carried  by  the  electrolyte.  At  the 
cathode-current  collector  the  conditions  remain  the  same.  These 
boundary  conditions  are  shown  in  columns  three  and  four  of 
Table  1. 


"  =  (1+-K1-2.)  O1  --)v,.r  +  2'.i:-gCs(1+,)),  (12) 
"r-oriofn^(“+’v“r“Jc‘n+’))'  (13) 

where  Ds  is  the  diffusion  coefficient  of  Li  in  the  particle;  R  is  the 
universal  gas  constant;  T  is  the  temperature;  aR,  or,  and  ah  are  the 
radial,  tangential,  and  hydrostatic  stresses;  E  and  v  are  the  Young’s 
modulus  and  Poisson’s  ratio  of  the  electrode  material.  The  Li  partial 
molar  volume,  Q,  accounts  for  the  swelling  of  the  particle  upon 
lithiation  and  the  stress-diffusion  coupling  increases  with  CL  The 
hydrostatic  stress  is  defined  as: 


ah  =  (aR  +  2aT)/3.  (14) 

At  the  particle  center,  r  =  0,  the  displacements  and  Li  flux  are 
zero,  ur  =  0,  J  =  0.  At  the  particle  surface,  r  =  Rs ,  Li  enters  the  particle 
and  mechanical  surface  pressure  is  applied. 


J  =  js  at  r  =  Rs, 

(15) 

ffR  =  Prmcro  at  r  =  Rs. 

(16) 

The  Li  pore  wall  flux  at  the  particle  surface,  js,  is  described  by  the 
Butler-Volmer  surface  kinetics  model  which  predicts  the  current / 
Li  flux  across  the  boundary  as  a  function  of  the  electric  potentials 
and  concentrations  in  both  phases: 
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js  =  lf\C 


s\r=Rs 1 


(cs,max 


“  cs|r=Rs 


I'O  =  ^2(q,max-Cl)“C(C1)“A, 

k<  •cw(v~u'(Ci "■-«•))• 


V  —  01  —  027 


,17x  homogenization  operators;  for  details,  the  reader  is  referred  to 
references  [15]  and  [2]. 

The  micro-scale  volumetric  strain,  ey  is  found  via: 


(18) 


'  Ur\r= 


r\r=Rs 


Ri 


) 


3 


(25) 


2.4.  Ragone  plot  simulation 


where  cs,max  and  ci>max  are  the  maximum  Li  and  Li+  concentrations 
in  the  electrode  particle  and  electrolyte;  a  a  and  ac  are  the  anodic 
and  cathodic  transfer  coefficients;  p  is  the  surface  overpotential; 
and  U'(cs]r=Rs)  is  the  open  circuit  potential  as  a  function  of  Li 
concentration. 

The  surface  pressure,  Pmicro»  has  contributions  due  to  the  in¬ 
clusion  of  the  particle  in  the  swelling  aggregate,  Ph,  and  from  sur¬ 
face  stress,  Pss,  which  can  be  written  for  spherical  particles  as  [24]: 

^rnicro  =  (19) 


Pss 


nTo  +  Ks  £r\r=Rs 

Rs 


(20) 


where  to  is  the  deformation-independent  surface  tension,  K*  is  the 
surface  modulus,  and  rj  is  the  tangential  strain,  which  in  spherical 
coordinates  is  defined  as:  ej  =  ur/r.  The  model  presented  here  does 
not  include  more  complex  phenomena  such  as  non-Fickian  diffu¬ 
sion,  the  effects  of  chemical  potential  on  stress,  finite  strains,  and 
plastic  deformations  [25-27]. 


The  set  of  equations  presented  in  Sections  2.1 -2.3  is  used  to 
predict  the  Ragone  curve  for  the  battery  by  simulating  the 
discharge  of  the  battery  over  a  range  of  rates.  We  discharge  the 
battery  using  a  LixC6  anode,  PE0-LiCF3S03  electrolyte,  LiyMn204 
cathode  system  where  0  <  x  <  1  and  0  <  y  <  1.  Material  properties 
are  given  in  Tables  2-4.  For  the  simulations  in  this  paper,  the 
separator  is  50  pm  thick  and  the  cathode  is  100  pm  thick.  This  yields 
a  theoretical  capacity  of  the  battery  of  34.8  Ah  m-2.  The  theoretical 
capacity  for  a  given  design  is  the  minimum  of  the  anode  and 
cathode  capacities  upon  discharge.  For  the  cathode,  the  theoretical 
capacity  is  the  amount  of  cathode  active  material  multiplied  by  the 
amount  of  Li  that  can  be  intercalated  into  the  cathode  particles: 

*  =  <5CC 

J  (1  —  £)dx(cs,max  —  Cs,o)*^-  (26) 

x  =  5sc 

Similarly  the  capacity  of  the  anode  is  the  amount  of  anode  active 
material  multiplied  by  the  amount  of  Li  present  in  the  anode  at  the 
start  of  the  discharge: 


2.3.  Governing  equations  at  meso-scale 


The  influence  of  the  micro-scale  particles  on  the  macro-scale  is 
seen  through  the  effective  pore  wall  flux,  jeff,  and  chemical  chem¬ 
ically  induced  eigenstrain  eCh-  Similarly,  the  microscopic  pore  wall 
flux,  js,  and  surface  pressure,  Ph  depend  on  both  macro-  and  micro¬ 
scale  state  variables.  Aggregate  theory  [28  and  Bruggeman  re¬ 
lations  [16,20]  are  used  to  relate  the  scales.  The  effective  pore  wall 
flux  is  related  to  the  micro-scale  pore  wall  flux  through: 


Jeff  = 


3(1-4, 

n  Js  • 


(21) 


The  porosity  of  the  electrode  influences  the  effective  properties 
Deff  and  xeff  which  are  found  via  Bruggeman  relations  based  on  the 
diffusion  coefficient  D\  and  electrolyte  conductivity,  /c^  in  the  bulk 
electrolyte: 


Deff  — 


Keff  =  £3/2K  00 


(22) 


The  micro-scale  surface  pressure  and  macro-scale  chemically 
induced  eigenstrain  depend  on  macro-  and  micro-scale  strains 
through  the  Mori-Tanaka  model  [28-30].  Flere,  we  assume  there  is 
no  chemically  induced  strain  in  the  electrolyte. 


Ph  =  (bs  +  BsCeff  (e  -  ech)  j  n, 


(23) 


eCh  =  (1  -  £)ev  +  (1  -  4  (cr1  -  Cf^bs,  (24) 

where  Q  and  Cs  are  the  isotropic  constitutive  tensors  of  the  elec¬ 
trolyte  and  active  particles  and  Ceff  is  the  effective  material  tensor  of 
the  composite  electrode.  The  matrices  Bs  and  bs  are 


x  =  <5  AS 

J  (1  -  £)dx*cSio*F.  (27) 

X  =  <5  AC 

The  thickness  of  the  anode  of  142  pm  was  chosen  so  that  both 
electrodes  would  have  the  same  theoretical  capacity.  The  battery  is 
discharged  at  a  constant  rate,  /,  until  the  potential  difference  be¬ 
tween  the  electrodes  (</>i  |x=rcc-0i  lx=rAC)  reaches  the  cutoff  po¬ 
tential,  0!  of  2.7  V;  this  cutoff  potential  was  selected  due  to  its 
position  on  the  discharge  curves  —  as  below  2.7  V,  the  battery 
potential  quickly  drops  to  zero  with  a  negligible  increase  in 
capacity. 

At  both  scales,  the  governing  equations  are  discretized  spatially 
using  finite  elements  and  in  time  using  an  implicit  Euler  backwards 
scheme.  An  adaptive  time  stepping  algorithm  is  used  that  varies  the 
time  step  size  based  on  the  rate  of  change  of  the  potential  differ¬ 
ence  across  the  electrodes  of  the  battery.  Spatially,  we  discretize  the 
system  using  20  elements  in  the  separator  region,  80  elements  in 
each  electrode  region,  and  20  elements  in  each  micro-scale  prob¬ 
lem.  Studies  were  performed  to  ensure  mesh  and  time 


Table  2 

Anode  material  parameters,  LixC6,  0  <  x  <  1. 


Symbol  Value  Units  Description 


Cs,max 

26,400 

mol  m  3 

Cs>  0 

13,070 

mol  m  3 

Vs 

100 

S  m1 

Ds 

5e-13 

m2  s"1 

Es 

15 

GPa 

Vs 

0.3 

- 

Q 

4.221e-6 

m3  mol  1 

k2 

le-10 

m4  mor1  s" 

Maximum  Li  concentration  in  cathode 
Initial  Li  concentration  in  cathode 
Conductivity  in  cathode 
Diffusivity  of  Li  in  cathode 
Young’s  modulus 
Poisson’s  ratio 
Partial  molar  volume 
Butler— Volmer  reaction  rate  constant 


Table  3 

PEO— UCF3SO3  electrolyte  material  parameters. 
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Symbol 

Value 

Units 

Description 

Cl, max 

3920 

mol  m  3 

Maximum  Li+  concentration  in  electrolyte 

Cl,0 

2070 

mol  m  3 

Initial  Li+  concentration  in  electrolyte 

Kao 

6.5e-3 

S  m1 

Conductivity  in  electrolyte 

Di 

7.5e-12 

m2  s"1 

Diffusivity  of  Li  in  electrolyte 

El 

0.69 

MPa 

Young’s  modulus 

V\ 

0.4 

- 

Poisson’s  ratio 

convergence.  Two-point  Gauss  integration  is  used  at  the  macro¬ 
scale  which  results  in  320  micro-scale  problems.  In  total,  there 
are  564  macro-scale  state  variables,  and  33,280  micro-scale  state 
variables. 

Fig.  3  is  a  Ragone  plot  of  the  rate  versus  usable  capacity  tradeoff 
seen  when  the  battery  is  discharged  at  14  rates  ranging  from 
0.01  A  m-2  to  25  A  m-2;  these  results  use  a  uniform  porosity  of 
e  =  0.3  and  particle  radii  of  5  pm  in  both  electrodes.  At  low 
discharge  rates,  the  battery  utilizes  up  to  99%  of  the  theoretical 
capacity  of  the  battery.  However  as  the  discharge  rate  increases, 
there  is  a  significant  drop  in  the  usable  capacity  of  the  battery. 
Ideally,  the  usable  capacity  of  the  battery  would  not  decrease  with 
increasing  discharge  rate. 

We  compare  the  discharge  of  the  battery  using  a  carbon  anode 
to  a  Li  foil  anode.  The  theoretical  capacity  of  the  battery  is  the  same, 
however  the  Li  foil  provides  an  effectively  unlimited  source  of  Li+  at 
the  anode-separator  boundary,  Tas-  As  seen  in  Fig.  3,  this  shifts  the 
Ragone  line  to  the  right,  i.e.  higher  capacities  for  higher  rates. 
Nevertheless,  the  decrease  in  usable  capacity  for  high  rates  is  still 
seen.  The  advantage  of  analyzing  the  capacity  of  the  half-cell  of  the 
cathode  cycled  vs.  a  Li  foil  anode  is  in  the  decrease  in  computational 
cost  as  the  anode  does  not  need  to  be  discretized.  For  the  compu¬ 
tational  model  used  in  this  study,  the  number  of  macro-scale  state 
variables  decreases  to  404  and  the  number  of  micro-scale  state 
variables  is  halved  from  33,280  to  16,640. 

When  discharging  the  half-cell,  the  cutoff  potential  is  set  at 
3.2  V.  The  higher  cutoff  potential  as  compared  to  the  carbon  anode 
simulations  is  due  to  the  constant  anode  potential  (versus 
increasing  potential)  during  discharge  for  the  half-cell  simulations. 
For  both  types  of  studies,  the  cutoff  potentials  were  selected  based 
on  reaching  the  regime  in  which  the  overall  battery  potential 
quickly  drops  to  zero.  Qualitatively,  the  Ragone  plots  predicted  by 
this  model  are  in  good  agreement  with  experimental  results  for  Li 
batteries  17]. 

3.  Design  optimization 

In  order  to  improve  the  high  rate  usable  capacity  of  the  battery, 
we  optimize  the  electrode  structure  when  the  battery  is  discharged 


Table  4 

Cathode  material  parameters,  LiyMn204,  0  <  y  <  1. 


Symbol 

Value 

Units 

Description 

Cs,max 

22,900 

mol  m  3 

Maximum  Li  concentration  in  cathode 

Cs.O 

4351 

mol  m  3 

Initial  Li  concentration  in  cathode 

ffs 

3.8 

S  m1 

Conductivity  in  cathode 

Ds 

7.08e-15 

m2  s-1 

Diffusivity  of  Li  in  cathode 

Es 

10 

GPa 

Young’s  modulus 

0.3 

- 

Poisson’s  ratio 

Q 

3.497e-6 

m3  mol-1 

Partial  molar  volume 

k2 

le-10 

m4  mol-1  s-1 

Butler— Volmer  reaction  rate  constant 

To 

-0.606 

J  m  3 

Surface  tension  (Si  value) 

K5 

10.65 

N  ITT1 

Surface  modulus  (Si  value) 

Fig.  3.  Ragone  Plot  for  full  battery  simulation  with  z  =  0.3  and  Rs  =  5  |im,  for  both 
types  of  anodes,  the  vertical  line  is  the  theoretical  capacity. 


over  a  range  of  discharge  rates,  Ijj  =  1  :Nj,  where  N,  is  the  number  of 
rates  simulated.  Previous  studies  have  shown  an  influence  of  the 
porosity  and  micro-scale  particle  radii  on  the  discharge  behavior  of 
the  battery  [15,31,2].  Therefore  we  consider  the  local  porosities  and 
electrode  particle  radii  in  both  electrodes  as  design  variables.  For  a 
given  layout  of  the  electrodes,  we  simulate  the  discharge  of  the 
battery  over  the  specified  range  to  the  cutoff  potential  and  maxi¬ 
mize  the  minimum  capacity  for  the  range  of  discharge  rates 
considered.  This  max-min  multi-objective  optimization  problem  is 
solved  via  a  bound  formulation  [32].  To  this  end,  we  introduce  an 
auxiliary  variable,  (3: 

P  =  min(/jf>),  (28) 

where  Tj  is  the  time  the  battery  reaches  the  cutoff  potential  when 
discharged  at  rate  Ij.  The  objective  is  to  maximize  /3  by  manipulating 
the  design  variables,  s.  Constraints  are  placed  on  the  difference, 
(3  -  IjTj,  in  order  to  ensure  that  /3  is  the  minimum  value.  This  results 
in  the  following  problem  formulation: 

max/?(s) 

s  '  ~  (29) 

s.t.Gj  =  (3  —  Ijtj  <  0  V  j  = 

The  design  variables  of  the  particle  radii  are  limited  to  be  be¬ 
tween  1  nm  and  10  pm  and  the  porosities  are  limited  to  between 
0.05  and  0.95  to  ensure  that  both  electrode  and  electrolyte  phases 
are  present  throughout  the  electrode.  In  addition,  we  place  con¬ 
straints  on  the  maximum  stress  level  in  the  micro-scale  electrode 
particles  during  discharge: 

G0  =  max  max  max<jTjtr  t  -  <jmax  <  0, 

Vj  =  1  :  Nj,  t  =  1  :Ntj,  i  =  1  :  Nm,  V  =  1  :  Nr  ’ 

where  r  is  the  node  id  of  the  particle  finite  element  model,  i  the 
individual  micro-scale  particle,  and  t  the  time  step.  Ntj  is  the 
number  of  time-steps  needed  when  discharged  at  rate  Ij,  Nm  is  the 
number  of  micro-scale  particles,  and  Nr  is  the  number  of  nodes  in 
the  particle  model.  In  this  paper,  we  set  amax  to  50  MPa.  This  stress 
level  was  chosen  on  the  basis  of  stress  levels  seen  in  the  constant 
porosity  solutions  of  Ref.  [2].  To  reduce  the  number  of  stress  con¬ 
straints,  we  approximate  the  maximum  particle  stress  in  an 
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Table  5 

GCMMA  parameters. 


Parameter 

Value 

Step  size,  As 

0.01-0.2 

Initial  adaptation  of  asymptotes 

0.5 

Adaptation  of  asymptotes 

0.7 

Maximum  number  of  sub-cycles 

1 

Relative  change  in  optimization  variables 

As  =  10“8 

individual  particle,  i,  at  discharge  rate  Ij  using  the  Kreisselmeier— 
Steinhauser  (KS)  function  with  oj  =  3: 

°Tj,i  =  7,,n5I  ^exP(MfTTj,t,r,i)  *  =  1  :  Ntj,r  =  0  :RS. 

0J  t  r 

(31) 

The  accuracy  of  this  approximation  has  been  monitored  in  the 
numerical  examples  to  ensure  the  approximation  is  accurate  to 
within  0.1%  at  50  MPa. 

The  design  optimization  problems  are  solved  using  the  Globally- 
Convergent  Method  of  Moving  Asymptotes  (GCMMA)  of  Svanberg 
[33  .  The  GCMMA  constructs  a  sequence  of  convex  separable  sub¬ 
problems  that  are  solved  by  a  primal-dual  method,  and  is  guaran¬ 
teed  to  converge  to  a  Karush-Kuhn-Tucker  (KKT)  optimal  point. 
This  algorithm  is  specifically  suited  for  problems  with  large 
numbers  of  design  variables  and  few  inequality  constraints  as  it  is 
the  case  in  the  current  study.  The  GCMMA  is  considered  converged 
if  the  constraints  are  satisfied  and  the  relative  change  in  design 
variables,  As  is  sufficiently  small: 


cm-l  I 


<  As 


(32) 


where  m  is  the  iteration  number  of  the  optimization  process.  Pa¬ 
rameters  used  by  the  GCMMA  are  given  in  Table  5;  for  a  detailed 
explanation  of  these  parameters,  the  reader  is  referred  to  Ref.  [33]. 

GCMMA  requires  the  derivatives  of  Gj  and  Ga,  with  respect  to  the 
design  variables.  The  derivatives  of  the  objective,  jS,  with  respect  to 
the  design  variables  are  zero,  while  the  derivatives  of  the  con¬ 
straints,  Gj  and  Ga,  are  nonzero  due  to  the  dependence  of  the 
simulation  time,  ?j,  and  model  behavior  on  the  design  variables. 
Due  to  the  potentially  high  number  of  design  variables  and 
computational  cost  of  the  model,  the  sensitivities  are  found  using 
an  adjoint  approach.  The  main  advantage  of  this  approach  is  that 
the  numerical  costs  of  computing  the  design  sensitivities  are  nearly 
independent  of  the  number  of  design  variables.  For  the  problem  of 
interest,  the  adjoint  sensitivity  equations  only  need  to  be  solved  NG 
times  in  each  step  of  the  optimization  procedure,  where  NG  is  the 
number  of  constraints.  This  involves  integrating  backward  in  time 
the  linearized  multi-scale  model  described  in  Section  2.  Details  on 
the  formulation,  numerical  procedure,  and  computational  effi¬ 
ciency  of  performing  adjoint  sensitivity  analysis  for  the  multi-scale 
framework  are  provided  in  Golmon  et  al.  [2]. 

4.  Design  optimization  examples 

4.1.  Optimization  of  cathode  discharged  against  Li  foil  anode 

We  solve  the  above  design  optimization  problem  (29)  for  a 
range  of  discharge  rates,  Ij  =  [0.01:20]  A  m-2,  first  with  a  Li  foil 
anode.  The  cathode  is  initialized  with  uniform  porosity  of  0.3  and 
particle  radii  of  5  pm,  consistent  with  the  nominal  design  in  pre¬ 
vious  studies  [15].  Fig.  4(a)  shows  the  effect  of  optimizing  the 
electrode  design  on  the  Ragone  plot,  and  Fig.  4(b)  and  (c)  plots  the 


Capacity  Ahm~2  distance  from  anode,  /im 

(a)  Ragone  Plot  (b)  Local  Porosity 


distance  from  anode,  pm 


(c)  Local  Particle  Radii,  pm 

Fig.  4.  Ragone  plot,  porosity  and  particle  radii  distributions  of  initial  and  optimal  designs  for  design  optimization  problem  with  Ij  =  [0.01:20]  A  m-2. 
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initial  and  optimal  porosity  and  particle  radii  distributions.  A  39% 
improvement  in  capacity  is  obtained  at  20  A  m~2,  however  at  low 
rates,  there  is  a  27%  decrease  in  capacity.  The  theoretical  capacity  of 
the  battery  decreases  from  34.80  Ah  m-2  to  29.8  Ah  m-2  due  to  the 
increased  porosity  of  the  electrode.  The  optimal  configuration 
achieves  75%  of  the  theoretical  capacity  at  20  A  rrr2  and  up  to  86% 
of  the  theoretical  capacity  at  lower  discharge  rates.  The  optimiza¬ 
tion  problem  is  dominated  by  the  performance  at  the  20  A  rrr2 
discharge  rate  which  has  the  lowest  capacity  and  highest  electrode 
particle  stress  levels  throughout  the  optimization  process.  This 
would  allow  for  a  reduction  of  the  optimization  problem  to  a  simple 
single-objective  problem,  maximizing  the  capacity  at  a  discharge 
rate  of  20  A  rrr2. 

The  influence  of  the  electrode  design  on  the  behavior  of  the 
battery  is  seen  in  Figs.  5  and  6  which  show  the  Li+  concentration  in 
the  electrolyte,  the  normalized  Li  concentration  in  the  electrode 


particles,  and  the  stress  levels  in  the  electrode  particle  as  a  function 
of  time  and  position  within  the  battery  for  initial  and  optimal 
configurations  for  10  A  nrr2  and  20  A  m-2  discharge  rates. 

Two  main  mechanisms  dominate  the  discharge  process:  (1) 
lithiation  of  the  cathode  particles  at  the  back  of  the  electrode  as 
seen  in  Fig.  5(d),  and  (2)  depletion  of  Li+  in  the  electrolyte  within 
the  cathode  as  seen  in  Figs.  5(a),  6(a)  and  6(b).  Depletion  of  Li+  in 
the  electrolyte  is  undesirable  as  it  prohibits  full  utilization  of  the 
active  material.  Once  the  electrolyte  is  depleted,  the  potential  of  the 
battery  drops  quickly.  The  optimal  design  for  these  rates  increases 
the  porosity  towards  the  back  of  the  electrode  which  allows  for 
increased  Li+  transport  in  the  electrolyte  and  thus  a  longer  time 
before  depletion  occurs.  The  electrode  particle  size  is  decreased 
near  the  separator-cathode  interface  primarily  to  decrease  the 
stress  levels  in  the  particles;  smaller  particles  allow  for  faster 
diffusion  of  Li  into  the  particles,  decreasing  the  stress  and  strain 


ci  (M) 


Q  (M) 


0  .  50  100  150 

distance  from  anode,  /im 

(a)  ci  (M),  Initial  Config 


distance  from  anode,  /im 
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Max.  gt  (MPa) 
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Fig.  5.  Q  and  cs  concentrations  and  maximum  aT  as  a  function  of  time  and  position  in  the  battery  for  initial  and  optimal  configurations  when  discharged  at  /  =  10  A  m  2  for  the 
optimization  problem  with  Ij  =  [0.01:20]  A  m-2.  The  vertical  lines  in  (a)  and  (b)  show  the  separator-cathode  interface. 
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Fig.  6.  Q  and  cs  concentrations  and  maximum  aT  values  as  a  function  of  time  and  position  in  the  battery  for  initial  and  optimal  configurations  when  discharged  at  /  =  20  A  m  2  for 
the  optimization  problem  with  Ij  =  [0.01:20]  A  m-2.  The  vertical  lines  in  (a)  and  (b)  show  the  separator-cathode  interface. 


mismatch  in  the  particle.  The  larger  particle  sizes  slow  the  rate  of 
Li+  uptake  by  the  particles,  which  results  in  enhanced  Li+  transport 
towards  the  back  of  the  cathode  [2  .  The  optimal  porosity  distri¬ 
bution  is  most  important  in  increasing  the  battery  capacity,  while 
adding  the  radii  as  design  variables  allows  for  a  further  improve¬ 
ment  on  the  order  of  0.5—1  Ah  rrr2. 

For  the  discharge  at  20  A  m-2  using  the  initial  configuration, 
only  the  front  half  of  the  cathode  is  saturated  with  Li  during 
discharge  (Fig.  6(c))  and  the  cathode  particles  at  the  current  col¬ 
lector  interface  only  see  a  negligible  increase  in  Li  concentration. 
The  stress  constraint  level  of  50  MPa  is  exceeded  near  the  sepa¬ 
rator-cathode  interface  (Fig.  6(e)). 

In  contrast,  for  the  optimal  configuration,  the  Li  concentration 
throughout  the  cathode  increases  during  discharge  (Fig.  6(d)) 
which  increases  the  capacity  and  percent  utilization  of  the  theo¬ 
retical  capacity  of  the  battery.  The  Li+  concentration  at  the  back  of 
the  cathode  (Fig.  6(b))  is  depleted  over  a  smaller  region  of  the 


cathode  and  occurs  later  in  the  discharge.  The  stress  levels  in  the 
particles  are  maintained  below  50  MPa  (Fig.  6(f))  with  the  highest 
stress  levels  positioned  at  the  largest  particle  sizes.  The  “waves” 
seen  in  the  contour  lines  for  the  normalized  Li  concentration  in  the 
cathode  particles  are  due  to  the  interplay  between  the  different  Li 
transport  rates  in  the  electrolyte  and  particles. 

If  the  range  of  discharge  rates  is  reduced  to  a  range  of 
[0.01 : 10]  A  m-2  or  [0.01 :1  ]  A  m-2,  the  capacity  of  the  battery  can  be 
increased  over  what  is  possible  for  the  [0.01 :20]  A  m-2  range.  Fig.  7 
shows  the  designs  and  Ragone  plots  for  all  three  optimization 
problems  as  compared  to  the  initial  configuration.  The  decrease  in 
the  rate  range  to  [0.01 : 10]  A  m-2  results  in  an  optimal  design  with  a 
theoretical  capacity  of  33.96  Ah  rrr2  and  the  [0.01:1]  A  rrr2  range 
optimal  design  has  a  theoretical  capacity  of  45.21  Ah  rrr2. 

At  low  discharge  rates,  the  low  optimal  porosities  result  in 
higher  usable  and  theoretical  capacity  of  the  battery;  here  the  ca¬ 
pacity  is  reaction-rate  limited.  As  the  discharge  rate  is  increased, 
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(a)  Ragone  Plot  (b)  Local  Porosity 


(c)  Local  Particle  Radii,  fim 


Fig.  7.  Ragone  plot,  porosity  and  particle  radii  distributions  of  initial  and  optimal  designs  for  all  three  design  optimization  problems. 


and  the  transport  is  limited  by  diffusion  of  Li+  in  the  electrolyte,  the 
porosity  is  increased  which  allows  for  greater  Li+  transport.  Particle 
size  primarily  influences  stress  levels  as  large  particles  are  diffusion 
limited  resulting  in  concentration  gradients  that  cause  mismatch 
strains  and  stresses  due  to  the  stress-diffusion  coupling;  see  Eqns. 
(12)  and  (13).  With  very  small  particles  (<0.5  pm),  the  stress  levels 
will  increase  due  to  the  surface  stress  effect;  see  Eq.  (20).  The 
particle  size  can  also  influence  the  usable  battery  capacity  by 
effecting  the  local  rate  of  Li+  uptake  throughout  the  electrode. 

A  tradeoff  is  seen  between  the  low  and  high  rate  performance  of 
the  battery:  for  the  highest  rates  an  increase  in  usable  capacity  is 
shown  while  the  low  rates  show  a  decrease  in  usable  capacity.  The 
highest  discharge  rate  dominates  the  design  for  each  simulation, 
suggesting  that  the  electrode  should  be  designed  based  on  the 
highest  discharge  rate  expected  during  operation.  Optimal  designs 
for  each  discharge  rate  give  upper  bounds  on  what  is  achievable 


Fig.  8.  Ragone  plot  for  the  cathode  found  via  optimization  of  the  cathode  with  a  Li  foil 
anode  when  discharged  with  a  carbon  anode. 


with  these  design  variables,  giving  the  Pareto-optimal  set  of  de¬ 
signs  (Fig.  7(a)). 

As  a  next  step,  we  replace  the  Li  foil  anode  with  a  carbon  anode 
as  in  Section  2.4,  which  has  the  optimal  cathode  porosity  and  radii 
distribution  from  the  Li  foil  [0.01 :20]  A  m-2  discharge  rate  range.  As 
seen  in  Fig.  8,  the  performance  improvement  is  significantly 
reduced  and  the  overall  theoretical  capacity  of  the  cell  is  still 
reduced  by  24%.  This  result  demonstrates  that  it  is  necessary  to 
perform  the  design  optimization  on  the  full  problem  if  the  carbon 
anode  is  to  replace  the  Li  foil  anode. 

4.2.  Optimization  including  carbon  anode 

The  design  optimization  problem  (29)  is  solved  for  discharge 
rates  of  Ij  =  [0.01:1]  A  rrr2,  [0.01:10]  A  nr2,  and  [0.01:20]  A  nr2, 
with  both  anode  and  cathode  porosities  and  radii  as  optimization 
variables.  Both  electrodes  are  initialized  with  uniform  porosities  of 
0.3  and  particle  radii  of  5  pm.  For  simplicity,  the  stress  constraint 
(30)  is  applied  only  in  the  cathode  region. 

Simultaneously  optimizing  the  porosity  and  particle  size  dis¬ 
tributions  in  a  full  battery  cell  results  in  greater  capacity 
improvement  than  using  the  optimal  cathode  design  with  the  Li  foil 
anode.  As  shown  in  Fig.  9,  the  resulting  optimal  design  for  the 
[0.01 :20]  A  nr  2  range  improves  the  capacity  of  the  battery  by  61% 
relative  to  the  initial  configuration  when  discharged  at  20  A  m  2, 
while  at  low  discharge  rates  the  capacity  is  decreased  by  42%. 
Similar  levels  of  improvement  are  seen  at  the  [0.01:1]  A  nrT2  and 
[0.01:10]  A  nr  2  rate  ranges.  The  resulting  optimal  designs  are 
depicted  in  Fig.  10.  This  figure  also  shows  the  optimal  cathode 
configuration  for  the  Li  foil  anode  optimization  problem  for  the 
[0.01:20]  A  m-2  discharge  rate  range. 

The  high  porosity  towards  the  separator  region  of  the  battery  in 
both  electrodes  is  due  to  the  limited  amount  of  Li  available  in  the 
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Fig.  9.  Ragone  plot  for  initial  and  optimal  configurations  for  the  carbon  anode  opti¬ 
mization  problems. 


and  electrolyte  in  the  cathode  than  the  anode,  resulting  in  a  larger 
spatial  variation  of  the  radii  in  the  cathode. 

To  further  understand  the  transport  limitations  of  the  battery,  in 
Fig.  11  we  examine  the  internal  Li+,  normalized  Li  concentration, 
and  stress  fields  within  the  battery  for  the  initial  and  optimal  de¬ 
signs  at  a  20  A  m-2  discharge  rate.  The  initial  configuration  limits 
the  reactions  in  the  electrodes  to  the  region  closest  to  the  separator; 
the  slow  Li+  transport  limits  the  amount  of  the  electrodes  utilized. 
The  optimal  configuration  shows  reactions  occurring  throughout 
both  electrodes,  although  Li+  depletion  is  still  the  limiting  factor.  In 
the  initial  configuration,  the  stress  levels  in  the  cathode  exceed  the 
stress  constraint  level  as  shown  in  Fig.  11(e).  The  stress  levels  in  the 
optimal  configuration  were  decreased  primarily  via  the  increased 
porosity  and  the  variation  in  particle  radii  enhances  Li+  transport 
through  the  depth  of  the  cathode. 

The  study  above  shows  that  optimizing  both  electrodes  simul¬ 
taneously  results  in  greater  improvements  in  usable  capacity  than 
optimizing  a  half  cell.  Note,  that  these  simulations  were  performed 
for  a  discharge  of  the  battery  and  result  in  asymmetric  designs 
between  the  anode  and  cathode,  optimization  of  a  full  charge/ 
discharge  cycle  could  change  the  optimal  configuration. 


anode  as  opposed  to  the  effectively  unlimited  Li  source  when  a  Li 
foil  anode  is  simulated.  Increasing  the  porosity  towards  the  sepa¬ 
rator  region  increases  the  Li+  transport  rates  in  the  electrolyte 
allowing  a  higher  percentage  of  the  electrodes  to  be  lithiated.  At  the 
highest  rate,  the  cathode  porosity  is  further  increased  at  the  back  of 
the  cathode  to  allow  for  increased  Li+  transport  in  the  electrolyte.  In 
the  anode,  the  optimal  particle  radii  distributions  result  in  a  small 
improvement  in  Li+  transport  in  the  electrolyte  while  in  the  cath¬ 
ode,  the  stress  constraint  also  influences  the  distribution.  There  is  a 
greater  difference  in  diffusion  coefficients  between  the  electrode 


(a)  Porosity 


(b)  Radii 

Fig.  10.  Initial  and  optimal  electrode  designs  for  dual  electrode  optimization.  The 
additional  dashed  lines  in  the  cathode  region  are  the  optimal  designs  for  the  Li  foil 
anode  optimization  problem. 


5.  Conclusion 

In  this  paper,  a  design  optimization  approach  based  on  a  multi¬ 
scale  battery  cell  model  has  been  presented  and  applied  to  opti¬ 
mizing  the  electrode  design,  in  order  to  increase  the  usable  capacity 
of  a  battery  over  a  range  of  discharge  rates.  Optimal  porosity  and 
radii  distributions  in  the  battery  electrodes  were  found  via  a  multi¬ 
objective  design  problem  with  constraints  placed  on  the  stress 
levels  in  the  cathode  particles.  Studies  were  performed  on  both  a  Li- 
foil  half  cell  and  carbon  anode  full  cell  configurations  which 
resulted  in  different  optimal  designs.  For  all  studies,  the  highest 
rate  discharge  dominated  the  performance  and  the  porosity  design 
variables  had  a  greater  influence  on  the  usable  capacity  than  the 
particle  radii  design  variables.  Simulations  presented  showed  that 
the  transport  of  Li+  in  the  electrolyte  was  the  dominate  process 
driving  the  optimal  designs  resulting  in  increasing  porosities  with 
increasing  rates.  The  particle  radii  design  variables  were  important 
in  stress-mitigation  near  the  separator-cathode  interface  and 
served  to  further  enhance  the  Li+  transport  due  to  the  effect  of 
particle  size  on  the  rate  at  which  Li+  enters  the  particles.  A  tradeoff 
was  observed  between  low  and  high  rates  because  the  highest  rate 
of  discharge  drives  the  optimization  problem. 

Design  optimization  of  the  electrode  porosities  and  particle  radii 
resulted  in  a  39%  increase  in  capacity  for  the  Li-foil  half  cell  and  a 
61%  increase  in  capacity  for  the  full  cell  versus  the  non-graded 
design  at  the  20  A  m  2  discharge  rate.  Flowever  the  optimization 
results  are  limited  by  the  level  of  accuracy  of  the  model  which  does 
not  include  effects  such  as  solid  electrolyte  interphase  (SEI)  layer 
growth,  chemical  side  reactions,  and  other  degradation  phenom¬ 
ena.  Only  a  single  discharge  of  the  battery  was  simulated,  and 
optimizing  the  electrode  configuration  over  a  full  charge/discharge 
cycle  may  result  in  different  optimal  electrode  structures.  Further 
optimization  studies  over  a  full  charge/discharge  cycle  and  over 
multiple  cycles,  including  degradation  mechanisms,  additional 
design  variables  of  the  electrode  thicknesses,  and  at  higher 
discharge  rates  are  recommended. 

Nevertheless,  the  results  presented  in  this  paper  suggest  that 
functionally  grading  the  electrode  can  improve  the  overall  battery 
performance  allowing  for  greater  usable  energy  density  in  a  bat¬ 
tery.  Furthermore,  this  study  has  demonstrated  the  utility  of 
computational  design  optimization  as  an  efficient  tool  for  the 
development  of  battery  cells.  As  our  optimization  approach  can 
incorporate  a  large  number  of  design  variables  and  can  be  extended 
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Fig.  11.  Q  and  cs  concentrations  and  maximum  Oj  as  a  function  of  time  and  position  in  the  battery  for  initial  and  optimal  configurations  when  discharged  at  /  =  20  A  m  2  for  the  full 
battery  optimization  problem.  The  vertical  lines  show  the  anode-separator  and  separator-cathode  interface. 


onto  more  complex  micro-  and  macro-scale  models,  this  method 
could  be  applied  to  simulations  investigating  the  overall  battery  cell 
design  such  as  those  in  Ref.  [11]. 
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